BIIL-MR-}9n 


I 


AD-A254  142 


MEMORANDUM  REPORT  BRL-MR-3992 


BRL 


AN  EXPERIMENTAL  FACILITY  FOR 
PACKED  AND  FLUIDIZED  BED  STUDIES 


CSABA  K.  ZOLTANI 


JULY  1992 


^UG101932|| 


APPROVED  FOR  PUBUC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


U.S.  ARMY  LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 


0  2  8 


92-21855 


I 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 

Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information  Service. 
U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield.  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position, 
unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report  does  not  constitute  indorsement 
of  any  commercial  product. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burd«n  lor  this  collection  of  information  is  estimated  to  average  i  hour  per  response,  including  the  time  for  reviewing  instruaions.  searching  ensting  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services.  Oireaorate  for  information  Operations  and  Reports.  12 IS  >efferson 
Oavts  Highway.  Suite  1204.  Arlington,  VA  22202-1302.  and  to  the  Office  of  Management  and  Rudget.  Paperwork  Reduction  Project  (0704-0188).  Washington,  DC  20S03- 


1.  AGENCY  USE  ONLY  (Leave  bUnk) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

July  1992 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final,  Jun  91  -  Feb  92 


5.  FUNDING  NUMBERS 


An  Experimental  Facility  for  Packed  and  Fluidized  Bed  Studies 


PR:  1L161102AH43 


6.  AUTHOR(S) 

C.  K.  Zoltani 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AOORESS(ES) 


B.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground.  MD  21005-5066 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


BSL-MR-3992 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

The  design  of  an  experimental  facility  for  the  study  of  drag  resistance  in  stationary  and  fluidized  gun 
propellant  beds  is  described.  It  establishes  the  capability  for  the  simulation  of  flows  at  five  times  the  Reynolds 
number  and  in  tubes  of  four  times  the  cross-sectional  area  heretofore  repented.  The  design  makes  possible  for 
the  first  time,  at  Reynolds  numbers  encountered  in  a  ballistic  packed  be^  unlimited  running  times.  In 
addition,  transient  flow  conditions  can  also  be  simulated. 


14.  SUBJEa  TERMS 


fluidized  bed  processes,  ballistics,  packed  beds,  flow  measurement,  Reynolds  number 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  _  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED  SAR 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev  2-89) 

Prescribed  by  ANSI  Std  ;38-i8 
298-102 


Intentionally  left  blank. 


TABLE  OF  CONTENTS 


Page 

LIST  OF  FIGURES  .  v 

LIST  OF  TABLES .  v 

ACKNOWLEDGMENTS  .  vii 

1.  INTRODUCTION  .  1 

2.  REQUIREMENTS  OF  THE  NEW  TEST  FACILITY .  3 

3.  DESIGN  OF  THE  FLOW  RIG  .  4 

3.1  Background .  4 

3.2  Specification  of  the  Components .  7 

3.2.1  Straight  Pipe  Section .  7 

3.2.2  Loss  at  an  Orifice .  11 

3.2.3  Flow  Straightener .  11 

3.2.4  Perforated  Plates  and  Baffles  .  13 

3.2.5  Screens  .  13 

3.2.6  Test  Bed  .  14 

3.2.7  Sonic  Nozzle  .  14 

3.2.8  The  Choice  for  a  High  Pressure  Gas  Source .  14 

3.2.9  Optical  Access  .  16 

3.2.10  Tube  Material  Strength  .  16 

4.  RESULTS  .  17 

5.  CONCLUSIONS  .  17 

6.  REFERENCES .  19 

APPENDIX:  THE  BLOWDOWN  PROCESS .  21 

LIST  OF  SYMBOLS  .  25 

DISTRIBUTION  LIST .  27 

iii 


Intentionally  left  blank. 


IV 


LIST  OF  FIGURES 


Figure  Page 

1.  Side  View  of  Flow  Rig  for  Fixed  Bed  Studies .  8 

2.  Instrumentation  for  Laser  Doppler  Anemometry  (LDA)  Velocity  Field  Analysis 

(Top  View) .  9 

3.  Particle  Image  Velocimetry  (PIV)  Diagnostics  for  Velocity  Measurements 

(Top  View) .  10 


LIST  OF  TABLES 


Table  Page 

1 .  Packed  Bed  Simulators .  1 

2.  Drag  Correlations  in  Use .  2 

3.  Pressure  Loss  Summary  .  15 


D7IC  QU  ALITY 


5 


I  Accession  For 

,  BTIS  GRAAI  ~ 

j  DTIC!  TAB 
j  Uri  '»r7  ir!'.,r)cf!d 


By - - - 

Dl  St-l’.'  hut, 

Ava  11,  bilit-’  'oies 


iDlst 


■’■'ll  aud/ur 
li eoial 


V 


Intentionally  left  blank. 


VI 


ACKNOWLEDGMENTS 


It  is  a  pleasure  to  thank  Mr.  Fred  Robbins,  Ballistic  Research  Laboratory  (BRL),  for  sharing  his 
experience  related  to  fixed  bed  resistance  experiments.  Thanks  are  also  due  to  Ms.  Liz  Marcou,  BRL,  for 
drawing  the  figures. 


Intentionally  left  blank. 


Vlll 


1.  INTRODUCTION 


Gun  perfonnance  predictions  rely  on  the  availability  of  packed  and  fluidized  bed  drag  data.  Due  to 
the  difficulty  and  cost,  only  a  very  limited  amount  of  exoerimental  bed  drag  data  is  available  at  the  present 
time.  To  fill  this  gap,  a  new  flow  rig,  which  allows  routine  data  generation  under  a  wide  variety  of 
conditions,  was  designed. 

The  resistance  to  flow  through  packed  beds  has  many  industrial  applications  and  has  received  wide 
attention.  Some  of  its  relevance  to  the  ballistic  problem  has  been  summarized  recently  in  Zoltani  (1991). 
Industrial  interest  usually  centers  on  low  velocity  flows  and,  therefore,  the  experimental  results  must  be 
interpreted  with  caution  when  used  in  a  ballistic  context.  Current  practice  has  its  antecedents  in  the 
research  of  Forchheimer  (1901)  and  Ergun  (1952)  and  is  basically  an  attempt  to  record  pressures  on  the 
walls  of  the  container  of  the  packed  bed.  Cold  flow  experiments  of  interest  to  baUisticians  have  been 
reported  by  Jones  and  Krier  (1983),  Wilcox  and  Krier  (1980),  Kuo  and  Nydegger  (1978)  and  Robbins  aixl 
Gough  (1978)  (see  Zoltani  [1991,  Tables  1  and  2]). 


Table  1.  Packed  Bed  Simulators 


Reference 

Bed 

Diameter 

(cm) 

Test  Section 
Length 
(cm) 

Particle 

Scale 

(cm) 

D/dp 

Pmu 

(MPa) 

Remark 

Ergun 

(1952) 

2.54 

20.0 

0.02-0.10 

25-111 

— 

glass,  lead 
shot 

Kuo  et  al. 
(1978) 

0.77 

30.0 

0.0826 

9.3 

14.0 

WC  870 
ball 

Robbins  et 
al.  (1978) 

7.62 

93.0 

1.39-30.3 

9.6-60 

20.0 

perf 

Jones  et  al. 
(1983) 

2.54-5.08 

20.3 

0.96-b.0 

8.5-50 

2.5 

glass  beads 

proposed 

work 

15.24 

100.0 

. . . 

1  cm  up 

15 

7.0 

— 

Table  2.  Drag  Correlations  in  Use 


Reference 

Coefficient  of  Drag,  Fv 

Shape 

Reynolds  No.A^oid  Regime 

Ergun 

(1952) 

F,  =  150  +  1.75  Ac 

various 

0.67  <  Re  <  2300 

0.4  <  e  <  0.65 

Kuo  & 

Nydegger 

(1978) 

F,  =  276.23  +  5.05(Re)°" 

spheres 

767  <  ^e  <  24  330 

0.37  <  e  <  0.39 

Robbins  & 

Gough 

(1978) 

F  =  2.5 

spheres 

cylinders 

Re  <  121  675 

0.396  <  e  <  0.509 

Jones  & 

Krier  (1983) 

F,  =  150  +  3.89(Re)°'" 

spheres 

737  <  Re  <  126  670 

0.38  <  e  <  0.44 

Note:  Here  fte  is  tne  Reynolds  number  divided  by  (1  -  e).  Note  that  the  correlation  of  Robbins  and  Gough  is  the  inverse  of 
F„  i.e.,  defined  as 

_  Ap  dp  gi 

where  the  symbols  are  the  same  as  defined  for  Equation  1.  For  details  see  Jones  and  Krier  (1983).  X  in  Robbins  and 
Gough  is  a  shqre  factor. 

In  these  experiments,  the  flow  was  assumed  to  be  steady  and  much  of  the  data  was  taken  at  moderate 
to  low  Reynolds  numbers.  The  flow  medium  being  bottled  gas,  the  time  available  for  the  acquisition  of 
data  was  limited.  This  fact  also  precluded  the  determination  of  velocity  profiles.  The  measurements  were 
used  to  establish  correlaticms  between  Reynolds  number  Re,  where  the  Reynolds  number  used  has  been 
divided  by  the  1  -  e,  and  friction  factor  for  beds  of  given  porosity,  see  Table  2.  There,  F^,  in  the 
formulation  of  Ergun  (1952),  is  the  friction  factor.  In  other  words, 


P  _  Ap  dp' 

“  L  pU  (i-e)^  ’ 


(1) 


where  Ap  is  the  pressure  drop,  dp  the  particle  diameter,  e  the  void  fraction  p  the  gas  viscosity,  and  U  the 

Re 

velocity  is  as  defined  in  Equation  3,  and  L  the  bed  length.  Values  of  F„  increase  with  increasing - 

1  -e 

from  150  into  the  thousands.  The  total  energy  loss  in  the  bed  is  the  sum  of  both  viscous  and  kinetic 
energy  losses,  but  depending  on  the  Re  number  regime,  one  or  the  other  may  predominate.  These 
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correlations  rely  to  various  degrees  on  Ergun’s  (1952)  results.  The  effect  of  the  proximity  of  the  walls 
is  not  directly  included.  Validity  of  the  correlations  are  restricted  to  the  indicated  Re  number  and  porosity 
regimes.  Extrsqwlation  of  the  formulas  to  outside  these  regimes  usually  entails  a  serious  error.  A  number 
of  unanswered  questions  remain,  including:  Are  radial  pressure  gradients  present  in  the  bed?  and  What 
is  the  nature  of  the  velocity  field  in  and  past  the  bed? 

The  next  sections  describe  the  requirements  on  the  experimental  facility  followed  by  the  actual 
specifications  required  to  attain  these  goals. 

2.  REQUIREMENTS  ON  THE  NEW  TEST  FACILITY 

To  extend  the  existing  database  on  packed  bed  propellant  drag,  an  experimental  rig  was  designed  to 
generate  flows  with  Reynolds  numbers,  based  on  propellant  length  scale  of  1  cm,  to  exceed  10^  in  the  test 
section.  Beds  up  to  155  mm  in  diameter  need  to  be  accommodated. 

The  experiments  must  yield  the  following  data: 

(1)  gas  mass  flow  rate 

(2)  pressure  at  several  stations  along  the  packed  bed 

(3)  temperature  at  the  entrance  and  exit  of  the  bed 

(4)  velocity  profiles  of  the  gas  flow  at  the  entrance  and  the  exit  of  the  bed 

(5)  boundary  layer  profiles  in  the  presence  of  ullage  and  in  the  two-phase  medium 

Propellants  of  different  shapes,  sizes  and  perforations  at  different  densities  of  loading  will  constitute  the 
bed  components.  Means  of  restraining  the  bed  from  moving  during  the  experiment  must  be  devised. 

The  measurement  of  velocity  profiles  in  tubes,  especially  near  curved  wall  surfaces,  is  a  challenging 
task.  For  an  adequate  resolution  of  the  boundary  layer  and  flow  profile  in  each  of  the  two  perpendicular 
coordinate  directions  at  each  cross  section  requires  between  15  and  20  data  points  per  radius.  The 
symmetry  of  the  flow  is  checked  by  taking  data  across  the  whole  cross  section.  The  data  taken  at  each 
and  every  location  must  be  repeatable  within  a  small  percentage  of  deviation.  Pradent  experimental 
practice  requires  that  data  be  checked  repeatedly  and  often  to  find  errors  in  the  flow  or  the 
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instrumentation.  The  bottom  line  is  that  at  the  upstrcam  and  downstream  locations  where  velocity  assays 
are  required,  in  excess  of  240  measurements  each  may  be  required  for  a  total  of  around  500  measurements. 

To  date,  propellant  bed  drag  experiments  were  run  with  bottled  gas.  However,  velocity  measurements 
impose  a  taxing  burden  on  the  gas  supply.  To  illustrate,  for  a  given  propellant  bed,  to  take  500  velocity 
measurements,  the  flow  taking  place  at  Re  =  10^  would  require  the  following  quantities  of  gas:  a  size  300 
gas  cylinder  of  1.75  ft^  capacity  (49.6  liters),  say,  initially  at  2,490  psi  (1.71  x  10’  Pa)  and  a  throat  area 
of  1.24  X  KT*  m^  (piping  of  5/32-inch  ID)  empties  in  about  174  seconds  when  the  flow  is  choked.  (For 
an  overview  of  the  blowdown  process  see  the  Appendix.)  In  a  reasonably  seeded  flow,  the  acquisition 
of  data  by  LDA  takes  about  10  seconds  per  point.  The  traversing  table  must  then  be  moved  to  the  new 
position,  which  takes  another  10  seconds  at  the  minimum.  Thus,  allowing  for  a  10%  wastage,  eight  data 
points  per  gas  cylinder  can  be  anticipated.  One  can  then  conclude  that  62  bottles  of  gas  at  $90  each  for 
ultra  high  purity  nitrogen  (since  it  must  be  oil  free),  or  about  $6,000  per  experiment  in  gas  cost  alone  can 
be  anticipated.  This  may  be  somewhat  optimistic  since  many  times,  due  to  uneven  seeding  and  other 
factors  related  to  the  nature  of  Laser  Doppler  Anemometry  (LDA)  data  acquisition,  more  time  per  data 
point  may  be  logged.  To  shut  off  the  gas  supply  between  data  positions  is  not  practical  since  flow 
transients  will  be  introduced  and  time  must  be  allotted  for  these  to  die  down.  One  sees  then  that  using 
bottled  gas  is  not  an  economical  way  to  proceed.  An  alternative  way  of  supplying  the  required  amounts 
of  gas  is  needed  and  is  suggested  below. 

3.  DESIGN  OF  THE  FLOW  RIG 

3.1  Background.  Drag,  or  resistance  to  flow,  can  be  inferred  from  the  pressure  drop  caused  by  the 
presence  of  the  obstacle  in  the  flow.  Differential  pressure  flow  metering  devices  such  as  orifices,  sonic 
nozzles,  and  venturis  rely  on  the  fact  that  the  pressure  distribution  in  a  flow  can  be  related  to  the  flow  rate. 
Thus,  if  the  pressure  drop  caused  by  the  presence  of  a  device  can  be  accurately  recorded,  a  direct 
relationship  to  the  flow  rate  can  be  established.  Accelerating  the  flow  through  a  constriction  of  known 
area  and  knowing  the  resulting  pressure  drop,  the  mass  flow  rate  can  be  determined. 

Three  coefficients  are  of  interest:  the  flow  coefficient,  Cj  ;  the  discharge  coefficient,  Cp  ;  and  the 
pressure  drop  coefficient,  K.  Recall  that  the  pressure  drop  between  two  axial  stations,  x,  and  Xj  in  a  tube 
under  incompressible  flow  assumptions,  is  directly  proportional  to  the  distance  Xj  -  x,. 
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(2) 


Ap/(0.5pU^)  =  f(x,  -  x,.^,Re) 


where  U  is  the  mean  approach  flow  velocity,  D  the  tube  diameter,  p  the  gas  density,  and  the  right  side 
is  a  functional  relationship.  We  note  that  the  Fanning  factor,  more  commonly  used  in  aerodynamics,  is 
one-quarter  as  large  as  f.  A  relationship  for  f  can  be  obtained  by  combining  the  continuity  equation 


m  =  pQ  =  pAU  , 


(3) 


with  Equation  2, 


1  _  p^U^A^ 
f  2pApA^ 


(4) 


Taking  the  square  root  of  both  sides  and  defining  f  = 


one  gets 


f  =  m/[A(2pApn  .  (5) 

Also,  it  should  be  noted  that  measurement  accuracy  is  increased  when  Ap  is  large.  Since  the  flow  rate 
Q  is  proportional  to  the  square  root  of  the  pressure  drop  (i.e.,  if  the  variation  in  p  is  10  to  1),  the 
corresponding  flow  rate  variation  is  3  to  1.  This  can  be  easily  inferred  from  m  =  pQ  =  fA^pAp  ,  i.e., 

Q  ~  t/^- 


For  a  flow  meter,  the  friction  factor  is  commonly  referred  to  as  a  flow  coefficient  Q  and  expressed 


as 


f  p 


^2Ap^ 


=  C  (X,  k/D,  Re)  , 


(6) 


5 


where  Q  is  the  volume  flow  rate.  X  is  the  area  ratio 


Ai,  k  is  the  roughness  of  the  pipe,  the  pipe  diameter  is  D, 


of  the  orifice  to  the  pipe  cross-sectional  area 
and  Re  is  the  Reynolds  number.  Note  that 


Q  =  f. 


Recalling  that  the  pressure  difference  upstream  of  the  constriction  to  that  at  the  minimum  cross- 
sectional  area  of  the  flow  can  be  related  to  the  area  ratio  as 


Ap-(1-X2).  (7) 

the  effect  of  the  magnitude  of  the  area  ratio  on  the  magnitude  of  Q  can  be  reduced  by  dividing  Ap  by 
(1  -  X*),  yielding  the  so-called  discharge  coefficient,  Cq. 

Co  =  (Q/A)/[p(l  -  X^)/(2Ap)]°'  .  (8) 

Co  is  the  true  flow  rate  divided  by  the  theoretical  flow  rate  and  the  relationship  between  Q  and  Co  is 
given  by 

Co  =  c,(l  -  .  (9) 

At  higher  Reynolds  numbers  the  discharge  coefficient  can  be  expressed  with  an  equation  of  the  form 

Cp  =  C_  +  b/Re "  ,  (10) 

where  the  C«  is  the  discharge  coefficient  at  infinite  Reynolds  number,  and  b  and  n  are  constants  (Miller 
1989).  Typically,  at  Re  =  10®,  the  discharge  coefficient  is  around  0.6.  Finally,  the  pressure  drop 
coefflcient  is  defined  as 

K  H  Ap/(0.5pU^)  =  [(!-(!)  X)m  ^)Y  .  (11) 

where  ^  is  the  coefficient  of  contraction  A^Aj;  A„  is  the  area  of  the  vena  contracta,  the  smallest  area  of 
the  discharge  past  the  orifice  plate;  Aj  is  the  orifice  cross-sectional  area;  and  X,  as  before,  is  an  area  ratio. 
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The  right-hand  side  of  Equation  1 1  follows  from  straightforward  manipulation  of  the  conservation  laws 
for  incompressible  flow. 

3.2  Specification  of  the  Components.  An  overview  of  the  flow  rig  is  given  in  Figure  1.  One  of  the 
innovations  introduced  is  that  the  rig  incorporates  its  own  gas  source,  furnishing  air  at  a  rate  high  enough 
to  ensure  that  the  desired  Reynolds  number  in  the  test  section  is  reached.  Upon  leaving  the  compressor, 
the  gas  enters  the  test  section  after  traversing  a  flow  straightener  (ASME-type,  perforated,  multiplate  flow 
rectifier).  The  gases  exiting  the  test  bed  go  through  a  second  flow  straightener  and  the  sonic  nozzle 
metering  section  before  exiting.  By  carefully  calibrating  the  sonic  nozzle,  it  may  be  adapted  to  function 
as  a  metering  device  both  in  the  choked  and  subsonic  regimes,  alleviating  the  need  for  an  orifice  plate  for 
flow  measurement  when  subsonic  conditions  prevail.  Pressure  probes  across  the  flow  meters,  along  the 
test  section  as  well  as  temperature  gages  fore  and  aft  of  the  test  bed  will  be  mounted  as  shown  in 
Figure  1.  Optical  access  to  the  flow  will  also  be  provided  to  enable  velocity  measurements  to  be  made 
by  state-of-the-art  flow  diagnostic  techniques  (Figures  2  and  3),  see  for  example  Dybbs  and  Ghorashi 
(1991). 

The  sonic  nozzle  flow  meter  assembly  of  Smith-Matz  geometry  with  a  1.798-inch  (0.045-m)  throat, 
contains  the  flow  rectifier,  inlet  temperature,  inlet  pressure,  and  throat  taps  with  a  balancing  manifold. 
Sonic  and  subsonic  calibration  need  to  be  provided  with  the  unchoking  test.  It  is  felt  that  the  discharge 
coefficient  variation  can  be  held  to  be  within  0.5%.  The  associated  instrument  package  includes  the  inlet 
temperature,  pressure  and  differential  pressure  probes,  and  transmitters.  The  data  acquisition  board  for 
a  286-type  computer  and  custom  software  will  display  the  mass,  volumetric  flow  rates,  and  the  status  and 
readings  of  the  pressure  probes. 

The  modus  operandi  adopted  here  was  to  start  at  the  back  of  the  rig  (exhaust  end),  calculate  the 
pressure  drop  through  each  of  the  components,  and  thus  arrive  at  the  performance  required  of  the  device 
furnishing  the  required  amount  of  air.  The  calculations  were  based  on  a  standard  6-inch  (152.4-mm)  pipe, 
which,  due  to  cost  considerations,  is  preferable  to  a  specially  fabricated  155-mm  tube. 

3.2.1  Straight  Pipe  Section.  The  pressure  drop  in  a  pipe  under  steady  flow  conditions,  as  already 
shown  in  Equation  2,  is  given  by  the  following  relationship, 

Ap  =  f(L/D)pUV2  ,  (12) 
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SCREEN  \  BAFFLE  PLATE 


Figure  1.  Side  View  of  Flow  Rig  for  Fixed  Bed  Studies. 
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Figure  2.  Instrumentation  for  Laser  Doppler  Anemometrv  (LDA)  Velocity  Field  Analysis  (Top  View). 
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Figure  3.  Particle  Image  Velocimetry  (PIV)  Diagnostics  for  Velocity  Mesurements  (Tod  View). 


where,  as  usual,  f  is  the  friction  factor,  L  the  length  of  the  pipe  section,  D  the  diameter,  and  p  the  density 
of  the  fluid.  Machining  costs  dictate  that  pipe  sections  should  not  exceed  1  meter  in  length. 
Commercially  available  steel  pipes  have  an  f  value  of  around  0.02  at  Re  =  10*.  Thus,  the  pressure  loss 
in  a  1 -meter  pipe  section,  under  the  conditiorts  of  interest  here  (p  =  1.2  kg/m*,  U  of  the  order  of  10  m/s, 
and  L/D  =  6)  will  be  185  Pa  and  the  K  value  is  0.12  based  on 

K  =  f  •  L/D  *  Ap/(0.5pU^)  .  (13) 

3.2.2  Loss  at  an  Oriflce.  For  a  Class  2  oriflce  (i.e.,  one  having  sharp  edges  and  made  of  thin  plate) 
and  t  <  O.ld  where  d  is  the  oriflce  opening  diameter,  t  the  oriflce  material  thickness,  and  Re  >  10*,  where 
the  Re  number  is  based  on  the  orifice  diameter,  in  the  nomenclature  of  Ward-Smith  (1980),  and  noting 
the  similarity  to  Equation  10,  i.e.,  ^  =  X, 


K  =  [(1- (t>^)/(4>^)]'  .  (14) 

Here  ^  is  the  orifice  opening  to  the  tube  area  ratio  and  (t>  the  contraction  ratio.  Ward-Smith  (1980)  gives 
tabulations  of  vs.  For  ^  =  0.5,  K  turns  out  to  be  4.954.  Equation  14  is  derived  by  using  the 
definitions  of  ^  in  conjunction  with  the  incompressible  conservation  equations  between  the  upstream, 
vena  contracta  and  downstream  positions  of  a  metering  device. 

3.2.3  Flow  Straightener.  The  accuracy  of  flow  metering  devices  depends  on  whether  the  flow  is  fully 
developed  and  on  the  success  of  eliminating  transverse  motion  (i.e.,  swirl)  in  the  flow.  Several  studies 
exist  on  the  effect  of  swirl  on  the  discharge  coefficient  of  an  orifice  plate  (Zanker  1962).  Kinghom  (1977) 
determined  discharge  coefficient  bias  errors  for  venturis  and  orifices.  For  a  swirl  angle  of  fifteen  degrees, 
typically,  the  error  can  approach  10%.  BSI  (1963)  stated,  as  was  also  confirmed  by  Kreith  and  Sonju 
(1965)  and  Miller  (1989),  that  the  persistence  of  swiil  may  extend  over  considerable  distances.  Without 
a  flow  straightener,  general  practice  suggests  (Miller  1989)  that  100  pipe  diameters  may  be  needed  to 
attenuate,  at  Re  =  10*  a  swirl  at  2°.  ISO  Standard  5167  (1980)  states  that  "Acceptable  velocity  profile 
conditions  can  be  presumed  to  prevail  when  at  each  point  across  the  pipe  cross  section  the  ratio  of  local 
axial  velocity  to  the  maximum  axial  velocity  at  the  cross  section  located  at  the  end  of  a  very  long  straight 
length  over  100  D  of  similar  pipe."  Flow  straighteners  remove  the  necessity  of  including  long  runs  of 
straight  tubing  to  achieve  fully  developed  flow  conditions. 
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Swill  has  several  effects  on  the  flow:  with  increasing  swiil,  the  position  of  the  maximum  axial 
velocity  moves  progressively  away  from  the  tube  axis.  Consequently,  there  is  a  reduction  of  the  axial 
velocity  along  the  axis  and.  at  high  enough  swiil,  reversed  flow  can  occur  at  the  pipe  axis.  Swiii  can 
defined  by  a  parameter  W  as  follows: 


W  =  4M_^/(ilpU^D^)  .  (15) 

where  is  the  angular  momentum  fiux.  As  shown  by  Baker  and  Sayre  (1974),  the  ratio  of  the  mean 
to  fully  developed  flow  without  swiil  friction  factors  can  be  given  as 

f„/f  =  1  +  11.5W  '  exp(0.45  L/D)  .  (16) 

where  L  is  the  distance  between  measuring  stations  and  D  the  pipe  diameter.  The  relation  is  valid  for 
1.25  X  10*  <  Re  <  2.0  x  10*.  Lugt  (1963)  describes  a  simple  device  which  can  indicate  the  presence  of 
swirl  in  a  pipe  flow.  It  consists  of  a  fitting  with  two  pressure  holes  side  by  side  on  the  periphery  of  the 
pipe.  A  swirl  at  the  wall  will  produce  a  differential  pressure  detectable  by  a  simple  manometer. 

Thus,  the  basic  objective  is  to  remove  major  flow  distortions  of  the  axial  velocity  so  that  fully 
developed  flow  conditions  exist  at  all  metering  devices.  Toward  these  ends,  a  number  of  authors  have 
suggested  remedies,  here  we  discuss  two  of  them.  It  is  important  to  minimize  the  value  of  K  of  these 
devices  so  that  the  pressure  loss  occurring  there  remains  small.  Sprenkle  and  Courtright  (1958)  proposed 
a  three-plate  arrangement  with  one  tube  diameter  spacing  between  plates  with  holes  amounting  to  about 
1/2  of  the  pipe  cross  section  distributed  over  the  plate.  A  value  of  K  =  16  was  measured  for  the  device. 
This  value  declined  to  1 1  when  the  inlet  side  of  the  orifices  was  beveled.  It  was  also  found  that  a  length 
of  eight  diameters  before  the  orifice  after  the  flow  straightener  assured  settling  of  the  flow  and  elimination 
of  metering  errors. 

An  improvement  in  the  K  value  was  reported  by  Zanker  (1962)  with  a  single  perforated  plate  but  with 
orifices  of  different  sizes  distributed  on  the  plate.  The  largest  hole  was  placed  near  the  axis  and  a 
honeycomb  (one  pipe  diameter  in  length)  was  attached  to  the  plate.  The  cells  were  square  with  their  axes 
coinciding  with  die  holes  in  the  plate.  The  settling  length  was  four  pipe  diameters  and  the  velocity  profile 
was  found  to  t^iproach  that  of  fully  developed  flow.  At  Re  =  10*.  the  K  value  was  5.7;  while  at  Re  =  10*, 
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it  declined  to  5.1.  Kinghom  et  al.  (1976)  report  a  further  improvement  with  a  straightener-venturi  nozzle 
package.  The  Mitsubishi  conditioner  (Miller  1989)  is  easier  to  fabricate  than  that  due  to  Zanker,  however, 
the  design  is  proprietary. 

3.2.4  Perforated  Plates  and  Baffles.  To  hold  the  test  bed  in  place,  either  a  screen  or  plate  and  a  baffle 
may  be  used.  A  perforated  plate  may  be  treated  as  an  orifice  plate  with  an  opening  equal  to  the  sum  of 
the  holes  in  the  plate. 

Assiune  a  plate  of  thickness  t  with  N  orifices  of  diameter  d.  The  porosity  of  the  plate  ^  is  defined 
as 

^s7tNdV(4A),  (17) 

where  A  is  the  pipe  cross  section.  For  t/d  <  0.8,  Ward-Smith  (1980)  gives  the  following  relation: 

K  =  [(0.609^(1  -  (t/d)^  (1  -  M+  ^^  *)-'  -  l]"  ,  (18) 

which  correlates  data  in  the  range  0.006  <  4  <  0-75  with  a  root-mean-square  error  of  10%  in  the  value  of 
K.  K  can  vary  widely;  indeed,  Ward-Smith  (1980)  quotes  values  from  0.57  to  3.58  x  10*. 

3.2.5  Screens.  The  upstream  end  of  the  propellant  bed  is  usually  held  in  place  by  a  screen.  A  wire 
mesh  may  be  characterized  by  the  number  of  wires  per  unit  length,  m  (units  of  inverse  length),  and  the 
porosity  of  wire  strands  of  diameter  d. 


^s(i-md)*.  (19) 

Gauzes  and  screens  offer  a  high  resistance  to  the  flow  compared  with  the  wall  shear  stress  in  developing 
or  fully  developed  flow.  A  sudden  change  in  pressure  is  experienced  when  flow  traverses  the  screen. 
Experimental  values  are  given  in  both  Ward-Smith  (1980)  and  Miller  (1990). 

Conventionally,  the  flow  through  each  aperture  in  a  screen  is  viewed  as  the  flow  through  a  single 
orifice  of  porosity  Note  that  Ward-Smith  (1980)  reports  that  at  hi^  Reynolds  numbers,  with  turbulent 
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flow,  K  is  independent  of  Reynolds  number.  Under  these  conditions,  de  Vahl-Davis  (1964)  gives  the 
following  relationship. 


K  =  K_+  55/Re  =  [( 1  -  +  55/Re  .  (20) 

K.  is  the  value  of  K  at  high  Reynolds.  Ward-Smith  (1980),  for  Re  >  10^  and  0.3  <  ^  <  0.9,  gives  curves 
of  K  as  a  function  of  a  pressure  parameter,  [(1  +  K_®*)^)‘‘.  Alternatively,  Ower  (1977)  gives  for  the 
resistance  of  a  gauze  screen  the  expression, 

K  =  6.5[(1  -  p)/p2][Ud/p(p/p)]->/'  .  (21) 

where  P  is  the  open  area  of  the  gauze,  U  the  velocity  ahead  of  the  screen,  d  the  diameter  of  the  wires,  and 
(p/p)  the  kinematic  viscosity  of  the  air.  No  general  formula  for  the  resistance  of  honeycombs  is  known, 
twt  some  data  is  available  from  Whitaker  et  al.  (1970)  and  Wieghart  (1953). 

3.2.6  Test  Bed.  The  pressure  drop  through  the  test  bed  was  calculated  on  the  basis  of  Ergun’s  (1952) 
theoiy  neglecting  the  first  term,  the  constant  150,  which  is  small  with  respect  to  the  retained  second  temi 
(see  Table  2). 


Ap/L  *  1.75[(1  -  e)/e^]pUVdp  .  (22) 

Here,  dp  is  the  particle  mean  diameter,  e  the  void  fraction,  and  p  the  gas  density. 

3.2.7  Sonic  Nozzle.  The  advantage  of  sonic  nozzles  as  flow  measuring  devices  is  that  upstream  of 
the  nozzle  changes  in  the  downstream  conditions,  as  long  as  the  nozzle  is  choked,  are  not  felt.  The 
pressure  recovery  across  the  nozzle  can  approach  92%,  however,  a  more  conservative  aixl  likely  figure 
is  85%.  It  is  also  possible  to  instrument  the  sonic  nozzle  so  that  it  can  act  as  an  orifice  plate  when  sonic 
condititms  are  not  reached. 

3.2.8  The  Choice  for  a  High  Pressure  Gas  Source.  The  estimated  pressure  drop  through  the  flow  rig 
is  given  in  Table  3.  The  Reynolds  number  at  the  entrance  to  the  test  section,  based  on  a  propellant 
particle  size  of  1  cm,  an  approach  velocity  of  6.41  m/s  and  gas  density  of  73.2  kg/m\  is  2.5  x  10^.  The 
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Table  3.  Pressure  Loss  Summary 


Component 

Pressure  Drop 
(Mpa) 

Conunents 

flow  straightener 

ASMS  type 

straight  section 

steel  tube 

wire  screen 

— 

test  section 

4.9340 

Ergun’s  eq. 

wire  screen 

0.0030 

— 

baffle  plate 

0.0315 

— 

flow  straightener 

0.0003 

ASME  type 

straight  section 

0.0003 

— 

sonic  nozzle 

0.2200 

est.  15%  loss 

particle  loading  e  was  taken  to  be  O.S.  Allowing  for  a  5%  pressure  loss  in  the  system  due  to  extraneous 
causes,  the  gas  must  be  supplied  at  a  pressure  exceeding  6.0  MPa. 

As  pointed  out  in  Section  2,  for  a  gas  source,  bottled  gas  in  view  of  the  anticipated  cost  is  not  a  viable 
option.  Rather,  an  situ  generation  of  the  pressurized  gas  suggests  itself.  Though  the  up-front 
investmentmay  be  higher  than  for  the  conventional  approach,  the  facility  will  pay  for  itself  in  gas  cost 
savings  in  a  short  period  of  time. 

Several  types  of  compressors  (see  Pichot  {1986]  and  Compressed  Air  and  Gas  Institute  [1989]  for 
details)  need  to  be  considered.  Basically,  one  can  distinguish  between  positive  displacement  types 
(i.e.,  reciprocating,  rotary,  dynamic,  centrifugal,  and  axial)  and  thermal  types,  that  is  injeaors. 
Reciprocating  compressors  can  produce  gases  in  the  tens  of  thousands  of  psig,  are  flexible  and  efficient. 
Their  Wg  drawback  is  that  they  require  large  foundations  due  to  the  unbalanced  inertial  forces  arvd,  most 
importantly,  the  gas  flow  is  discontinuous.  Thus,  their  use  in  the  packed  bed  studies  can  only  be 
comtemplated  in  conjunction  with  the  use  of  a  gas  accumulator. 

Screw-type  rotary  compressors  arc  capable  of  handling  flows  ranging  from  200  to  20,0(X)  m^/hr  under 
discharge  pressure  ranging  from  3  bar  gage  in  a  single  stage  to  about  13  bar  gage  in  two  stages.  These 
compressors  use  two  spiral  lobed  rotors  whose  hollows  are  filled  with  gas  at  the  inlet  port.  Compression 
begins  because  of  the  spiral  configuration  (i.e.,  during  rotation  the  contact  surfaces  move  parallel  to  the 


15 


axes  of  the  rotors  toward  the  outlet  causing  gradual  compression  until  the  trapped  air  is  put  at  the  outlet 
port). 

Centrifugal  compressors  can  reach  discharge  pressures  of  10,000  psi  (68.94  MPa)  with  an  inlet 
capacity  better  than  200,000  acfin  (94.38  mVs).  Staging  is  usually  employed.  Ambient  air  is  accelerated 
by  the  first  impeller,  then  a  diffuser  converts  the  velocity  into  pressure.  A  heat  exchanger  cools  the  air 
before  it  flows  through  a  second  stage  impeller  and  is  then  discharged  through  an  aftercooler  to  the  user. 
Axial  machines,  on  the  other  hand,  have  even  higher  capacities  but  their  discharge  is  limited  to  500  psi 
(3.44  MPa). 

3.2.9  Optical  Access.  Velocity  measurement  in  the  flow  requires  optical  access.  At  elevated 
pressures,  windows  must  be  specially  fabricated.  With  a  Young’s  modulus  of  50  x  10*.  modulus  of 
rupture  up  to  1(X),000  psi  (689.4  MPa),  sapphire  has  the  needed  strength  and  optical  qualities.  At  0.5  tun 
wavelength,  transmittance  approaches  85%  and  the  index  of  refraction  is  1.77.  The  hydrostatic  burst 
strength  of  simply  supported  windows  is  8,300  psi  (57.22  MPa).  Thus,  to  maintain  1,000  psi  (6.89  MPa) 
with  a  safety  factor  of  four,  the  manufacturer.  Union  Carbide  (1988),  recommends  a  thickness  to  diameter 
ratio  of  a  window  of  0.15. 

3.2.10  Tube  Material  Strength.  In  tubes  where  the  radius  to  the  tube  wall  thickness  exceeds  10,  the 
circumferential  stress  can  be  calculated  based  on  the  following  relation  (Gere  and  Timoshenko  1990), 

a,  =  (pr)/5  ,  (23) 

where  p  is  the  pressure  inside  tlie  tube,  r  the  radius,  and  t  the  wall  thickness.  The  longitudinal  stress 
follows  from 

Oj  =  (pr)/(2t)  .  (24) 

Taking,  as  a  typical  example,  a  steel  tube  of  0.15  m  (6  inch)  ID,  wall  thickness  of  0.0063  m  (0.25  inch), 
with  a  safety  factor  of  4,  say  at  13.78  MPa  (2,000  psi)  in  the  flow,  the  stresses  are  well  below  the 
allowable  yield  and  tensile  stresses  of  commercially  available,  resistance-welded  carbon  steel  pipes  of  68 
and  76  ksi  (4.68  x  10*  Pa),  respectively. 
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4.  RESULTS 


To  create  the  desired  flow  conditions  in  the  test  cell,  the  best  solution  is  to  use  a  compressor. 
Experiments  can  be  run  indefinitely,  a  great  advantage  when  performing  difficult  flow  surveys.  Unsteady 
flow  conditions  can  also  be  simulated. 

In  choosing  a  compressor,  several  trade-offs  must  be  considered.  Most  of  the  centrifugal  compressors 
on  the  market  guarantee  high  flow  rates  but  at  moderate  to  low  pressures.  The  present  requirement  calls 
for  high  pressures  by  compressor  standards,  but  at  fairly  low  flow  rates.  Also,  as  the  compressor  gas  exit 
velocity  increases,  due  to  the  U^  term  (Equation  22),  the  pressure  drop  through  the  bed  increases  rapidly 
putting  a  large  demand  on  the  gas  generator. 

To  reach  the  required  Reynolds  number  at  the  entrance  to  the  test  section,  and  assuming  gas  exit 
temperatures  of  300  K  (which  is  realizable  if  the  compressor  comes  with  a  heat  exchanger),  a  multi-stage 
compressor  is  required.  For  example,  the  Bauer  Model  B  28.3  three-stage  compressor  discharges  air  at 
3(X)  K,  a  density  of  73.20  kg/m^  and  a  pressure  of  6.40  x  10*  Pa.  At  a  flow  rate  of  250  cfm  (0.1 17  mVs), 
the  velocity  is  6.41  m/s,  yielding  a  Reynolds  number  of  2.5  x  10*. 

The  facility  will  allow  the  generation  or  data  heretofore  unavailable  and  needed  by  interior  ballistic 
modelers.  By  including  a  particle  disperser  and  a  means  of  recovering  the  propellant  in  the  experimental 
rig,  many  aspects  of  the  fluidized  bed  can  also  be  simulated.  This  could  include  but  is  not  limited  to 
ullage,  center  core,  and  velocity  field  data  in  addition  to  detailed  drag  data. 

5.  CONCLUSIONS 

A  new  experimental  facility  for  the  measurement  of  inert  propellant  bed  drag  has  been  designed.  Its 
operational  envelope  encompasses  Reynolds  numbers  in  excess  of  2.5  x  10*  in  test  beds  up  to  155  mm 
in  diameter.  Both  steady-state  and  time-dependent  operations  are  possible.  In  steady  state  there  is  no 
restriction  on  the  time  for  the  experiment  This  allows  the  taking  of  velocity  profiles  of  the  flow.  The 
experiment  is  also  instrumented  with  both  LDA  and  PIV  (particle  image  velocimetry)  data  acquisition 
systems  enabling  for  the  first  time  the  determination  of  both  the  pressure  along  the  test  section  as  well  as 
velocity  profiles  and  turbulence  data  fore  and  aft  of  the  packed  bed. 
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APPENDIX: 

THE  BLOWDOWN  PROCESS 
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Intentionally  left  blank. 
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THE  BLOWDOWN  PROCESS 


A  mass  M  of  gas  is  contained  in  a  tank  of  volume  V.  Through  an  opening  of  area  A  the  gas  exits. 
Choked  conditions  prevail  at  the  exit  from  time  t  =  0  to  time  t  =  tj.  It  is  desired  to  calculate  the  pressure 
history  during  the  emptying  process  as  well  as  the  time  needed  from  the  initial  opening  of  the  valve  to 
the  time  that  choking  conditions  can  no  longer  be  maintained.  It  is  assumed  that  the  flow  process  is 
isentropic  and  k  is  the  ratio  of  specific  heats  in  the  gas. 
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where  i  denotes  the  initial  state.  Differentiating  the  above  expression. 
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where  the  critical  mass  flux  is  given  by 
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Note  that  the  conditions  in  the  reservoir  are  continuously  changing  and  these  time-dependent  conditions 
have  to  be  reflected  in  the  use  of  Equation  (A-4). 
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Substituting  and  rearranging. 
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this  is  a  differential  equation  with  the  initial  conditions. 

at  t  =  0  ,  p/Pi  =  1  . 

Let  (p^  =  *  and  let  n  =  (1  -  3k)/2k.  The  equation  then  can  be  written  in  the  simplified  form  as  follows: 

X  "  =  C  .  (A-6) 

dt 

Here  C  is  a  constant,  the  value  being  the  right-hand  side  of  Equation  A-S.  This  is  the  classical  Bernoulli 
equation.  The  substitution  v  =  x"  *  ‘  reduces  it  to  a  linear  equation  which  was  first  solved  by  Leitmiz  in 
16%.  The  result  is 
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which  alternatively  can  be  solved  for  t  for  a  given  (p/p,)  ratio. 
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LIST  OF  SYMBOLS 


A  cross-sectional  area 
Cd  discharge  coefficient 
Q  flow  coefficient 
D  bed  (pipe)  diameter 
d  wire  diameter,  orifice  diameter 
dp  panicle  diameter 
F  drag  force 

f  friction  factor 

G  mass  flow  rate  per  unit  area 

g  acceleration  due  to  gravity 

h  hydraulic  height 

K  pressure  drop  coefficient 
k  roughness,  also  ratio  of  specific  heats 
L  length 

M  angular  momentum  flux,  also  total  mass 
m  number  of  wires  per  unit  length,  also  pQ 
N  number  of  particles 
p  pressure 

Q  volume  flow  rate 

Re  Reynolds  number  =  P^— 
t  plate  thickness 

U  mean  approach  flow  velocity 

V  volume 

W  swirl  parameter 

GREEK 

p  open  area 

e  void  fraction 

X  area  ratio 

^  viscosity 
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4  porosity  of  (date 
p  density 
o  stress 
X  shear  stress 
^  coefficiem  of  contraction 


SUBSCRIPT 


ang  angular 

c  critical  condition  (choked) 
i  initial  value 

m  mean 

p  [Kuticle 

V  velocity 

1  circumferential 

2  orifice 

oo  free  stream  condition 
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University  Park,  PA  16802-7501 
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1  Pennsylvania  State  University 
Assistant  Professor 

Department  of  Mechanical  Engineering 
ATTN:  Dr.  Stefan  T.  Thynell 
219  Hallowell  Building 
UnivCTsity  Park,  PA  16802-7501 

1  Pennsylvania  State  University 

Director,  Gas  Dynamics  Laboratory 
Department  of  Mechanical  Engineering 
ATTN:  Dr.  Gary  S.  Settles 
303  Mechanical  Engineering  Building 
University  Park,  PA  16802-7501 

1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Technical  Library 
333  Ravenwood  Ave. 

Menlo  Park,  CA  94025-3493 

1  Rensselaer  Polytechnic  Institute 
Department  of  Mathematics 
Troy,  NY  12I8I 

2  Director 

Los  Alamos  Scientific  Laboratory 
ATTN:  T3,  D.  Butler 

M.  Division,  B.  Craig 
P.O.  Box  1663 
Los  Alamos,  NM  87544 

1  General  Applied  Sciences  Laboratory 

ATTN:  J.  Erdos 
77  Raynor  Ave. 

Ronkonkama,  NY  1 1779-6649 

1  Battelle  PNL 

ATTN:  Mr.  Mark  Gamich 
P.O.  Box  999 
Richland,  WA  99352 

1  Stevens  Institute  of  Technology 

Davidson  Laboratory 
ATTN:  R.  McAlevy  III 
Castle  Point  Station 
Hoboken,  NJ  07030-5907 


1  Rutgers  University 

Department  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  S.  Temkin 
University  Heights  Campus 
New  Brunswick,  NJ  08903 

1  University  of  Southern  California 

Mechanical  Engineering  Department 
ATTN:  0HE200,  M.  Gerstein 
Los  Angeles,  CA  90089-5199 

1  University  of  Utah 

Department  of  Chemical  Engineering 

ATTN:  A.  Baer 

Salt  Lake  City,  UT  84112-1194 

1  Washington  State  University 

Department  of  Mechanical  Engineering 
ATTN:  C.  T.  Crowe 
Pullman,  WA  99163-5201 

1  Alliant  Techsystems,  Inc. 

ATTN:  R.  E.  Tompkins 
MN38-3300 
5700  Smetana  Drive 
Minnetonka,  MN  55343 

1  Alliant  Techsystems,  Inc. 

ATTN:  J.  Kennedy 
7225  Nmhland  Drive 
Brooklyn  Park,  MN  55428 

1  Science  Applications,  Inc. 

ATTN:  R.  B.  Edelman 
23146  Cumorah  Crest  Drive 
Woodland  Hills.  CA  91364-3710 

1  Battelle  Columbus  Laboratories 

ATTN:  Mr.  Victor  Levin 
505  King  Ave. 

Columbus.  OH  43201-2693 

1  Allegheny  Ballistics  Laboratory 

Propulsion  Technology  Department 
Hercules  Aerospace  Company 
ATTN:  Mr.  Thomas  F.  Farabaugh 
P.O.  Box  210 
Rocket  Center,  WV  26726 
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1  MBR  Research  Inc. 

ATTN:  Dr.  Moshe  Ben-Reuven 
601  Ewing  Sl,  Suite  C-22 
Princeton,  NJ  08540 


Aberdeen  Proving  Ground 


1  Cdr,  USACSTA 

ATTN:  STECS-PO,  R.  Hendricksen 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  uixlertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-MR-3992 _ Date  of  Report  July  1992 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct  Address 
in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DO-T 

Aberdeen  Proving  Ground,  MD  21005-5066 


OFRCIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

RRST  CLASS  ffRMT  Nb  0001,  M>6,  MO 


Postage  will  be  paid  by  addressee 


Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-5066 


